arXiv:l509.02219v2 [hep-ph] 10 May 2016 


Prepared for submission to JHEP 


AMS-02 fits Dark Matter 


Csaba Balazs’^ Tong LP 

ARC Centre of Excellence for Particle Physics at the Tera-scale, School of Physics and Astron¬ 
omy, Monash University, Melbourne, Victoria 3800 Australia 

E-mail: csaba.balazs@monash.edu, tong.li@monash.edu 

Abstract: In this work we perform a eomprehensive statistieal analysis of the AMS-02 
eleetron, positron fluxes and the antiproton-to-proton ratio in the eontext of a simplified 
dark matter model. We inelude known, standard astrophysieal sourees and a dark matter 
eomponent in the eosmie ray injeetion speetra. To prediet the AMS-02 observables we 
use propagation parameters extraeted from observed fluxes of heavier nuelei and the low 
energy part of the AMS-02 data. We assume that the dark matter partiele is a Majorana 
fermion eoupling to third generation fermions via a spin-0 mediator, and annihilating to 
multiple ehannels at onee. The simultaneous presenee of various annihilation ehannels 
provides the dark matter model with additional flexibility, and this enables us to simul¬ 
taneously fit all eosmie ray speetra using a simple partiele physies model and eoherent 
astrophysieal assumptions. Our results indieate that AMS-02 observations are not only 
eonsistent with the dark matter hypothesis within the uneertainties, but adding a dark mat¬ 
ter eontribution improves the fit to the data. Assuming, however, that dark matter is solely 
responsible for this improvement of the fit, it is diffieult to evade the latest CMB limits in 
this model. 
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1 Introduction 


Charged cosmic rays carry a wealth of information about galactic astrophysics and pos¬ 
sibly about new fundamental particle physics. Deciphering this information is, however, 
challenging because it requires the detailed understanding the injection and propagation 
of cosmic rays within the Galaxy. Fortunately, the last decade witnessed an increasing 
precision both in the experimental determination and the theoretical prediction of cosmic 
ray fluxes. As observations became more and more precise a deviation between them and 
prediction became apparent in the electron and positron fluxes [1-16]. The latest and most 
precise measurements of the electron, positron flux, antiproton-to-proton ratio, and proton 
flux came from the AMS-02 collaboration [17-21]. The increase of the positron spectral 
index and the growth of the positron fraction above 100 GeV are unexpected features of 
these measurements [17, 18]. 

The difference between these measurements and various predictions is the subject of 
debate. It may originate from unsatisfactory understanding of cosmic ray propagation, 
through unaccounted standard astrophysical sources (such as pulsars and/or supernova 
remnants), to more exotic new physics (such as dark matter annihilation) [22-24]. Moti¬ 
vated by the exciting possibility that the apparent excess of cosmic electrons and positrons 
is due to dark matter annihilation, in this work we examine whether the AMS-02 data are 
consistent with a typical particle dark matter model. First, we make a prediction for the 
expected background based on the propagation parameters of heavier cosmic isotopes and 
commonly used injection spectra. Then we calculate the contribution of dark matter an¬ 
nihilation to the electron, positron and anti-proton fluxes. Adding this to the background 
flux allows us to constrain the parameter space of the dark matter model. 

To determine the cosmic ray background due to standard astrophysical sources we 
adopt the following strategy. We assume that the relevant cosmic ray propagation pa¬ 
rameters and injection spectra can be determined by fitting the observed fluxes and the 
secondary-to-primary ratios of heavier nuclei (e.g. B/C,^° Be/®Be) and the low energy 
regions of the and p/p spectra. Based on these fits we derive the background for the 
and p/p fluxes. Then we calculate the injection spectra of and p due to dark mat¬ 
ter annihilation. Using the earlier determined diffusion parameters we propagate the dark 
matter annihilation products through the Galaxy. This procedure ensures a consistent as¬ 
trophysical treatment of cosmic rays originating from standard astrophysical sources and 
from dark matter. 

As particle physics description of dark matter we use the simplified model framework. 
This ansatz uses minimal and general theoretical assumptions. We consider a single dark 
matter particle, a Majorana fermion, that couples to standard fermions via a spin-0 media¬ 
tor. We do not assume a specific, single annihilation final state for the dark matter particle. 
Rather, more realistically and in line with minimal flavor violation [25], we allow the dark 
matter particle to annihilate into the third generation quarks and the tau lepton. The si¬ 
multaneous presence of various annihilation channels provides the dark matter model with 
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considerable flexibility, whieh enables us to simultaneously fit all oosmie ray speetra using 
a single partiele physies model and eoherent astrophysieal assumptions. This is one of the 
most important results of our work. Beyond this outcome we also delineate the AMS-02 
preferred region in the parameter spaee of the dark matter model. 

This paper is organized as follows. In See. 2 we describe the propagation equation 
and injection spectra for cosmic ray in galaxy. The values of eorresponding parameters 
are also given. In See. 3, we briefly deseribe the simplified dark matter model we use. Our 
numerieal results are given in See. 4. Finally, in See. 5 we summarize our main results. 


2 Injection and Propagation of Cosmic Rays 


Cosmie rays are energetie partieles propagating within the Galaxy, and are divided into 
primary and seeondary types [26-29]. Primary eosmie rays are likely to originate from 
powerful astrophysieal proeesses, sueh as supernova explosions and pulsars. By interaet- 
ing with intergalaetie matter they ereate seeondary oosmie rays [7, 8, 27, 29-31]. Propa¬ 
gation of oharged oosmie rays within the Galaxy oan be quantified by the diffusion model 
[32-35]. This model provides a meohanism to explain the retention and isotropio distribu¬ 
tion of high energy charged particles within the Galaxy, by desoribing particle scattering 
on Galactic media, sueh as magnetio fields [26, 29, 35, 36]. The speotrum of oosmie rays 
is modified by various energy loss meohanisms (due to interaction with the interstellar 
medium) and re-aeoeleration (due to interstellar shooks) [35, 37, 38]. 

Cosmie ray propagation within the galaetio halo is desoribed by the transport equation 
[35] 
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Here t,p) is the density of cosmie rays per unit of total particle momentum p, V is 
the convection velocity, and Tf{Tr) is the time soale for fragmentation (radioactive decay). 
The spatial diffusion coefficient is written in the form 


Dxx — f^Do{R/R q)^ , 


( 2 . 2 ) 


with R and (3 being the rigidity and partiele veloeity divided by light speed respectively. 
The diffusion coefficient in momentum spaee, i.e. Dpp, is proportional to the square of 
the Alfven veloeity va- The height of the eylindrieal diffusion halo is zq. The above key 
propagation parameters ean be eonstrained by fitting the seeondary-to-primary ratios of 
nuelei, that is the Boron-to-Carbon ratio (B/C) and the Beryllium ratio (^°Be/®Be). We 
adopt the diffusion re-aeoeleration model and the values of propagation parameters shown 
in Table 1, determined by the B/C and ^°Be/®Be data [39]. 
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Each cosmic ray species is deseribed by an equation as Eq. (2.1), with speeies speeifie 
parameters. The souree term of eosmie ray speeies i ean be generally deseribed by the 
produet of spatial distribution and injeetion speetrum funetions 

Qi{r,p) = f{r,z)qi{p). (2.3) 

Eor the spatial distribution of the injeeted primary eosmie rays we use the following su¬ 
pernova remnants distribution 

where the distanee between the Sun and the Galactie eenter is r© = 8.5 kpc, the height of 
the Galaetie disk is Zg = 0.2 kpc, and the two parameters a and b are taken to be 1.25 and 
3.56, respeetively. The normalization parameter /q is determined by the EGRET gamma 
ray data [40]. We assume the following power law with one break for the injeetion speetra 
of various nuelei 


Qi oc 


iR/Rir\R<Rl 

{R/Rir ^, R > Rl 


(2.5) 


and two breaks for primary eleetrons, i.e. R^^i, i?br2 with 71,72,73 being the power law 
indexs. 

Eollowing the approaeh in Ref. [39] we adopt a seale faetor Ce+ = 3.1 to take into 
aeeount the uneertainty in the ealeulation of the seeondary fluxes from proton-proton eol- 
lision eross section and enhancement faetor from heavier nuelei. It is introdueed to reseale 
the ealeulated seeondary flux to fit the data. The eorresponding injeetion parameters ean be 
determined by fitting the AMS-02 proton, eleetron, and positron data. We adopt injeetion 
parameters obtained by sueh a fit in Ref. [39]. The values of these injeetion parameters are 
shown in Table 1 . 

We use the Eisk potential {i = e~,e'^,p,p), relating the loeal interstellar fluxes 
to the one measured at the top of the atmosphere, to aeeount for the solar modulation 
effeet. We treat as speeies speeifie nuisanee parameters. Their best fit values are shown 
in Table 1. Sinee solar modulation affects the observed fluxes only below 10 GeV, the 
values of these parameters have no effeet on our eonclusions drawn about the dark matter 
eontribution. 


3 The Dark Matter Model 

In this seetion, we deseribe the partiele physies model we use to demonstrate that the 
AMS-02 data ean be explained by dark matter annihilation. In the reeent literature it 
was shown that Majorana fermions are one of the most plausible dark matter eandidates 
[41-48]. Inspired by this, we assume that dark matter is eomposed of Majorana fermion 
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propagation 

value 


nucleon injection 

value 


electron injection 

value 


solar modulation 

value 

Do (10^* cnF 

6.58 



1.811 


7i 

1.463 


A- (MV) 

1550 

<5 

0.33 


V2 

2.402 


72 

2.977 


(MV) 

1800 

Ro (GV) 

4 


Rl (GV) 

12.88 


73 

2.604 


0P (MV) 

518 

va (km 

37.8 


Ap (see caption) 

4.613 


(GV) 

2.858 


h (MV) 

0 

zo (kpc) 

4.7 


- 

- 


Rir2 (GV) 

68.865 


- 

- 

- 

- 


- 

- 


Ag (see caption) 

1.585 


- 

- 


Table 1. Parameters of propagation, nucleon/electron injection and solar modulation and their 
values adopted in our numerical analysis. The proton (electron) flux is normalized to Ap (AA at 
100 (25) GeV in the units of 10“® cm“^ sr“^ MeV“^. 


particles, which we denote by x- Motivated by the Higgs portal mechanism, we assume 
that the dark matter particle couples to standard fermions via a spin-0 mediator, that we 
denote by S [49, 50]. We cast the dark matter to mediator coupling in the form 

^ ^-^XlbXS. (3.1) 

Coupling between the dark matter and mediator is fixed to = 1. (This choice effectively 
absorbs into the mediator-standard model couplings.) Coupling between the mediator 
and standard model fermions / is given by 


CsDXfffS. (3.2) 

We assume that S only couples to third generation fermions, consistently with minimal 
flavor violation, i.e. / = b,t,T [25]. For simplicity we do not consider dark matter 
annihilation into a pair of S particles. With the interactions defined by Eqs. (3.1) and 
(3.2) dark matter annihilation is not velocity suppressed [51]. At the same time the dark 
matter-nucleon elastic scattering cross section is spin-independent (SI) and momentum 
suppressed. 

Under the above assumptions the dark matter model is described by the following 
parameters: 


P = {m^,ms,Xb,Xt,XT} ■ (3.3) 

The scan ranges for these parameters are 

1 TeV < < 10 TeV, 1 GeV < ms < 1 TeV, 10“^ < Xb, Xr, A* < 10^.(3.4) 

The potentially large values of the above effective couplings can only be understood in 
an underlying theory. They may include the effect of large but renormalizable perturba¬ 
tive couplings, large loop contributions from vector-like matter, resonant or Sommerfeld 
enhancements, or the combination of more than one such a factor [52]. 
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The source term arising from dark matter annihilation contributing to the cosmic ray 
species i is given by 


Q}{r,p) 


dNl \ 

2m2 ^ dE j ' 


(3.5) 


where {av) is the velocity averaged dark matter annihilation cross section, Bf = (av) // (av) 
is the annihilation fraction into the // final state, and diV/ /dE is the energy spectrum of 
cosmic ray particle i produced in the annihilation channel into //. In the parenthesis on 
the right hand side the total differential yield is the Bf weighted sum of the partial dif¬ 
ferential yields into specific final states. The sum includes contributions from all the third 
generation charged fermions ( 6 , f, r). AMS-02 plays an important role in constraining the 
coupling of the mediator to these fermions since Bj directly depends on these couplings. 

We use a generalized Navarro-Frenk-White (NFW) profile to describe dark matter 
spatial distribution within the Galaxy [53] 


Px(^) 


{r/rs) 

^\l + r/rsf-E 


(3.6) 


Here the normalization coefficient is po = 0.26 GeV/cm^ and the radius of the galactic 
diffusion disk is = 20 kpc. We fix the inner slope of the halo profile to 7 = 1. 


4 Results 


As discussed in Sec. 2, the propagation and injection parameters of cosmic rays are de¬ 
termined by fitting the B/C and ^°Be/®Be data and recent charged cosmic ray data from 
AMS-02, respectively [39]. The parameters in Table 1 thus imply prediction for cosmic 
ray measurements inferred from standard astrophysical sources. One can investigate the 
constraint on extra sources, such as dark matter, based on this fiducial astrophysical back¬ 
ground. 

To this end the Lagrangian of the dark matter model described in the previous section 
was coded in FeynRules [54]. Using model files generated by FeynRules, the annihilation 
fraction Bf and differential yields diV //dE in Eq. (3.5) were calculated by a modified 
version of micrOmegas 3.6.9 [55]. These dark matter model dependent variables were then 
input into the public code Galprop v54 [34, 37, 56-58] to ensure that near Earth cosmic 
ray fluxes from dark matter annihilation and background spectra obtained in a consistent 
way. 

The calculated cosmic ray fluxes, together with the measured spectral data points, 
were entered in a composite likelihood function, defined as 


21 n£ = ^ 


- fr? 

a? 


(4.1) 
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Here are the theoretieal predietions and are the eorresponding eentral value of 
the experimental data. The uneertainty cxj eombines the theoretieal and experimental un- 
eertainties in quadrature. We stipulate a 50% uneertainty of the theoretieal predietion 
of electron flux, positron flux and antiproton-proton ratio according to the estimates of 
Refs. [59-62]. This uncertainty takes into account, amongst other, the uncertainty related 
to the fixed propagation parameters. The sum in Eq. (4.1) runs over all the AMS cos¬ 
mic ray spectral data points: the electron flux (73 points), positron flux (72 points) and 
antiproton-proton ratio (30 points). We do not include the AMS-02 positron fraction data 
in the likelihood function; consequently the theoretical positron fraction flux is a prediction 
in our framework. 

Including observables from dark matter abundance, direct detection, or collider pro¬ 
duction in the likelihood function would not change its value significantly. We found that 
in the parameter region that dark matter annihilation can appreciably contribute to the 
charged cosmic ray fluxes the self-annihilation rate is high enough to decrease dark matter 
abundance below the observed level. In this case, assuming that x is just a component of 
dark matter, the likelihood is not affected by abundance. Dark matter direct detection is 
impaired by momentum suppressed x-nucleon elastic scattering cross section and the very 
high mass of x- As for the Large Hadron Collider (LHC), in the relevant parameter re¬ 
gion X particles are too heavy to produce in significant numbers via 14 TeV proton-proton 
collisions. 

LIG. 1 shows our main results: AMS-02 cosmic ray flux observations are consistent 
with the dark matter hypothesis within the uncertainties. The four frames display the 
various cosmic ray fluxes AMS-02 observed: electron flux, positron flux, positron fraction, 
and antiproton-to-proton ratio. AMS-02 central value measurements are shown by red 
dots and dark error bars indicate their uncertainty. The green solid line, on each frame, is 
obtained using the parameters shown in Table 1 and displays the predicted background flux 
originating from standard astrophysical sources. The blue solid line shows the prediction 
of the total cosmic ray flux with dark matter parameter values that best fit the AMS-02 
data. The blue curve is the sum of the background flux (green curve) and the dark matter 
contribution at the best fit point (magenta curve). A series of orange colored dots (forming 
vertical bars) indicate the theoretical uncertainty of the dark matter prediction given by the 
95% confidence region of dark matter model parameters. 

As the plots show adding a dark matter contribution to the background flux yields a 
better fit to the AMS-02 data. As expected, the electron flux is hardly changed by the dark 
matter contribution, while the latter somewhat improves the agreement between the theo¬ 
retical prediction and the antiproton-to-proton ratio data. This indicates that the dark mat¬ 
ter model is consistent with these data. The fit to the positron data is noticeably improved 
that justifies the addition of the dark matter component. Our likelihood function used to 
extract the best-fit dark matter parameters does not include the positron fraction data, that 
is the dark matter model parameters are not fit to the / (e"*" -f e~) fraction. Rather, after 
we extract the best fit dark matter model parameters, we calculate the positron fraction 
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Figure 1. Electron flux, positron flux, positron fraction, and antiproton-to-proton ratio observed by 
AMS-02 (red dots and dark error bars). The blue solid line shows the prediction of the total cosmic 
ray flux with dark matter parameter values that best fit the AMS-02 data. The total predicted flux 
is the sum of the background flux (green solid line) and the dark matter contribution. Orange dots 
indicate the 95% confidence region of the prediction. The magenta line is the flux from dark matter 
at the best fit point. 

using the best fit parameters. As shown by the blue eurve the e"*"/(e+ + e~) fraction data 
and the best fit (obtained without this data) agree very well. This is an important cross 
check of the internal consistency of the dark matter model and our parameter extraction 
procedure. 

The top frames of FIG. 2 show the regions of the dark matter parameter space pre¬ 
ferred by the AMS-02 data. Solid circles and squares denote the estimated 68% and 95% 
confidence regions, respectively. The favored mass of the dark matter particle is heavier 
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Figure 2. The AMS-02 favored region of masses (top left, ms vs. m^), couplings (top right, 
Xt/Xb vs. Xr/Xb), and cross sections (bottom, av vs. m^) in the simplified dark matter model we 
consider. The solid circles and squares estimate 68% and 95% confidence regions, respecfively. 
The besf fif poinf is indicafed by a friangle. 

than 2 TeV (at about 68% C.L.) with best fit point indicating an 9.3 TeV dark matter mass. 
The AMS-02 data favor a spin-0 mediator mass in the region of 1-700 GeV (at about 68% 

C.L.). 

For the mediator-SM fermion couplings the favored region indicates that the tau lepton 
coupling At- is generally larger than quark couplings Xb, Xt, being 1000 (10) times larger 
than Afe (A*) at the best fit point. This trend is governed by the electron and positron data fit: 
dark matter annihilations should produce mostly leptons to explain the difference between 
the astrophysical background and the AMS-02 data at high energies. The antiproton-to- 
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proton ratio data, on the other hand, require the moderate presenee of either bottom or top 
quarks in the final state. Henee the diagonal shape of the estimated 68% and 95% C.L. 
regions on the right hand frame of FIG. 2. The best fit point favors coupling values for 
which At ~ lOAi ~ lOOOAb. 

The bottom frame of FIG. 2 shows that the AMS-02 data require an effective dark 
matter annihilation cross section in the region of 1 x 10-23 _ 2 X 10-22 (5 X 10-24 - 
3 X 10-2^) cm^/s at about 68 (95) % C.L. An effective cross section so much higher 
than the standard thermal rate could indicate the non-thermal origin of self-annihilating 
dark matter particles responsible for AMS-02 [63-65]. Alternatively, the positron ray flux 
might receive a boost from dark matter substructure, such as over dense clumps, clouds, 
or disks which would allow for a reduced annihilation rate [66-75]. 

According to Ref. [76] a 1-10 TeV dark matter particle with an annihilation cross 
section of av ~ 10-2^ — IO -22 cm^/s, and dominant final state of or bb, is excluded 
by Planck and by Fermi-LAT gamma ray bounds from dwarf satellite galaxies. Since the 
annihilation rate at the recombination time places a (particle physics) model independent 
limit on the present day annihilation rate, either of these limits are hard to evade. Som- 
merfeld enhancement does not alleviate the problem, since the average relative velocity 
of scattering dark matter particles at the time of CMB is lower than the present day one. 
Uncertainties in the relevant astrophysical measurements, such as in the power injected 
into the CMB or the Fermi-LAT statistical/systematic errors, do not seem to leave enough 
room for the high dark matter annihilation cross section required to account for AMS-02. 
The most straightforward way to evade the Planck and Fermi-LAT limits appears to be 
including a standard, but presently unanticipated, astrophysical contribution to explain the 
AMS-02 measurements. With such additional contribution the dark matter annihilation 
cross section can be lowered and the model be made consistent with all data. 


5 Conclusions 

In this work we examined the plausibility of dark matter annihilation contributing to the 
recent AMS-02 data, the electron, positron fluxes and antiproton-to-proton ratio. On the 
top of the standard astrophysical cosmic ray flux prediction we included a dark matter 
component. Our choice of the dark matter model was a Majorana fermion coupling to third 
generation fermions via a spin-0 mediator. The initial flux from standard astrophysical 
sources and dark matter annihilation were propagated through the Galaxy using the same 
set of diffusion parameters. The latter were determined by fitting the cosmic ray fluxes of 
heavier elements and the low energy regions of the AMS-02 data. 

We have shown that not only AMS-02 observations are consistent with the dark mat¬ 
ter hypothesis within the uncertainties, but adding a dark matter contribution to the back¬ 
ground flux yields a better fit to the data. We also estimated the most plausible parameter 
regions of the dark matter parameter space in light of AMS-02. The observations prefer 
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a dark matter (mediator) mass in the 2-10 TeV (1-700 GeV) region at about 68% eonfi- 
denee level. The data also favor a dominant tau lepton-dark matter eoupling A,-, about ten 
times larger than top quark-dark matter eoupling \t at the best fit point. The antiproton- 
to-proton ratio data require that dark matter annihilation to quarks is dominated by either 
the top or the bottom final state with a slight preference for the latter. 

At the meantime we found it to be difficult to evade the CMB and Fermi-LAT gamma 
ray limits in this model due to the high annihilation cross section. With additional contri¬ 
bution to the positron spectrum from standard, but presently unknown, astrophysics this 
cross section can be lowered and the model be made consistent with all data. 
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